Mon. Not. R. Astron. Soc. 000, 000-000 () Printed 31 August 2009 (MN stylo file v2.2) 



Late-time supernova light curves: The effect of internal 
conversion and Auger electrons 
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ABSTRACT 

Energy release from radioactive decays contributes significantly to supernova light 
curves. Previous works, which considered the energy deposited by 7-rays and positrons 
produced by ^^Ni, ^^Co,^^Ni, ^''Co, ^^Ti and ''^Sc, have been quite successful in ex- 
plaining the light curves of both core collapse and thermonuclear supernovae. We 
point out that Auger and internal conversion electrons together with the associated 
X-ray cascade, constitute an additional heat source. When a supernova is transparent 
to 7-rays, these electrons can contribute significantly to light curves for reasonable 
nucleosynthetic yields. In particular, the electrons emitted in the decay of ^^Co, which 
are largely due to internal conversion from a fortuitously low- lying 3/2~ state in 
the daughter ^"^Fe, constitute an additional significant energy deposition channel. We 
show that when the heating by these electrons is accounted for, a slow-down in the 
lightcurve of SN 1998bw is naturally obtained for typical hypernova nucleosynthetic 
yields. Additionally, we show that for generic Type la supernova yields, the Auger 
electrons emitted in the ground-state to ground-state electron capture decay of ^^Fe 
exceed the energy released by the ^4^^ decay chain for many years after the explosion. 
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1 INTRODUCTION 



One of the great successes of nuclear astrophysics was 
demonstrating that the light curves of both thermonuclear 
(Type la) and, at least for late times, core collapse (Type 
Ib/c, Type II) supernovae are powered by the decay of ra- 
dionuclides synthesized in the explosion. Quickly following 
the dynamic explosion and nucleosynthesis phase a super- 
nova enters a state of homologous expansion. The temper- 
ature drops during the expansion so that nuclear fusion 
ceases and radioa ctive decay is t h e only means of cha nging 
the composition. PankevI l|l962t) , iTruran et al.l (Il967l ) and 
IColgate fc McKee ( 19691 ) first proposed the now widely ac- 
cepted view that the energy liberated in the decay of ra- 
dioactive ^''Ni followed by the decay of ^®Co to stable ^""Fe 
is the most important source for re-heating the supernova 
ejecta to temperatures high enough to shine brightly in the 
optical part of the electromagnetic spectrum. Initially, the 
bulk of the heating is induced by the energetic 7-rays pro- 
duced in the decays which thermalize and deposit their en- 
ergy via a number of processes, including Compton scat- 
tering and photoelectric absorption. Due to the expansion, 
however, the column density decreases with time as t~^, and 



the ejecta becomes more and more transparent to these high 
energy photons. 

^®Co decays predominantly via electron capture, but 
19% of all decays proceed via positron emission. The funda- 
mental importa n ce of this positron channel was first pointed 
out by lArnetj (|l979l '). Initially it was assumed that the 
positrons deposit all of their kinetic ener gy before form - 
ing positronium and annihilating locally l|Axelrodl Il980l ) . 
With that assumption, for a generic Type la supernovae 
(SN la) model, the ejecta are sufficiently transparent to 
7-rays after about ~ 200 days that the ^^Co positrons 
dominate the energy depo sition rate (e.g. IWooslev et al.l 
1 19891 : ISoUerman et al.ir2004l ) . Subsequently, the idea of com- 
plete and local ene r gy d eposition has been questioned. 
IChan fc Lingenfelteij ( 19931 ) suggested that as early as a few 
hundred days after the explosion positrons may begin to 
escape, especially i n the presence of radially combed mag- 
netic fields (see also Milne et al.l 199i ). However, it has also 
been argued (e.g. ISwartz et al.lll995 ), that the fraction of 
positrons which escape cannot be very large since the ob- 
served light curves could then not be explained without an 
additional source of heating. 

SN 1987A has demonstrated that light curves can be 
constructed for nearby supernovae even several years af- 
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ter the explosion. This has led to the consideration of less 
abundant but longer lived radioactive species. The most 
important of these is ^^Co, which is expected to be pro- 
duced in significant amounts as the decay product of the 
short lived ^'^Ni. Due to its relatively longer half-life (271.79 
days) and the higher opacity of the ejecta to the emit- 
ted 7-rays, it may dominate the bolometric light curve, 
especially o f core coll apse supe rnovac, at late times (e.g. 
IPinto et al.l[r988 : Wooslev et al.lil989. 'l. It is widely argued, 
that the only other radionuclide which n oticeably con- 
tributes to the bolometric light curve is **Ti (|Wooslev et al.l 
ll980l : lKumagai et al.lll989l'l. '^''Ti has an even longer half-life 
of 58.9 years ( Ahmad et al ]l2006'), and the short lived daugh- 



ter Sc has a strong positron channel, which means that the 
'''*Ti decay chain can be the dominant energy source a t very 
late times (Kuma gai et al.lll989l : IWooslev et al.lll989l ). 

In section [2] we review the physics of radioactive decay, 
emphasizing the roles of internal conversion (IC) and Auger 
electron production. We point out that, while ^'^Co has no 
positron channel, due to the nuclear structure of the daugh- 
ter ^^Fe, a significant amount of energy is released in the 
form of IC and Auger electrons. In section |3] we summa- 
rize the current observational status of late time supernova 
bolometric light curves of SNe la and SN 1998bw, a well 
observed member of the broad-lined SNe Ic class. For both, 
we sketch the effect electrons emitted in radioactive decays 
may have. Implications of our study are discussed and our 
conclusions are drawn in Section |4l 



2 RADIOACTIVE DECAY 

Most of the mass synthesized by fusion processes in super- 
nova explosions consists of nuclear species in the mass range 
A K 12 — 70. The nuclei are typically either stable or lie 
on the proton-rich side of stability. For the relevant nu- 
clei, fission and a-decay plays no role. The radioactive decay 
considered here occurs along an isobar towards more neu- 
tron rich nuclides in one of two ways, electron capture or 
positron emission. The first proceeds via the capture of an 
inner atomic (typically K or L shell) electron by a nuclear 
proton and the emission of an electron neutrino. The latter 
proceeds by the decay of a nuclear proton into a neutron 
and the emission of a positron and an electron neutrino. For 
both positron emission and electron capture, the transition 
to the daughter is a statistical process to a distribution of 
(excited) nuclear levels, with branching ratios given by the 
transition probabilities. 

Following electron capture, the daughter will have an 
electron hole in its atomic structure. Electrons cascade down 
to fill the gaps in the lower lying atomic shells and typically a 
series of characteristic fluorescence X-rays is emitted in the 
process. However, at every atomic transition, there is the 
possibility that instead of emission of a photon, the excess 
energy is carried away by one or several ejected outer elec- 
trons (Auger electrons). For ground-state to ground-state 
electron capture transitions, such as in the decay of ^^Fe, 
no 7-rays or positrons are emitted. Assuming that the neu- 
trino escapes without interactions, in such transitions the 
Auger electrons and the X-rays constitute the only sources 
of radioactive heating. 

An excited state of the daughter nucleus typically tran- 



sitions to lower lying states by emitting 7-ray photons until 
the ground state is reached. For each transition, however, 
there exists the additional possibility to de-excite via IC. 
In this process, the energy difference of the nuclear levels 
is carried away by the ejection of an inner atomic electron 
and there is no 7-ray photon emitted. These electrons are 
called IC electrons. The probability to de-excite via IC, as 
measured by the IC coefficient a = de-excitations 

^ #or 7 de — excitations ' 

generally small, but increases for low lying states. In the 
decay of ^^Co in particular, IC electrons are relatively copi- 
ously produced, due to the existence of a low-lying nuclear 
level in the daughter '^''Fe. For this 14.4 keV 3/2" state of 
^^Fe, IC is in fact the preferred mode of de-excitation, with 
Q = 8.58. 



3 LATE-TIME BOLOMETRIC LIGHT CURVES 

Late-time bolometric light curves have already been recon- 
struc ted from mult i -band photometry for severa l objects 
(e.g. ISuntzefj |2002| : ISoUerman et~al] |2002| . l2004h . For re- 
liable reconstruction, the contribution of the UV/optical 
iUBVRI) and near-infrared (JHK) bands should be in- 
cluded (UVOIR light curve). However, given that near- 
IR observations are rare at very late epochs, sometimes 
only _B-through-7 band observations are used, applying a 
near-IR correction ex trapolated from earlier epochs (e.g. 
ISoUerman et al.ll2002l ). 

Apart from the heating by various radioactive nuclei, a 
number of additional effects impact the slope of late-time SN 
light curves. Most importa ntly the ejecta become inc reas- 
ingly transparent to 7-rays jLeibundgut fc Pintol 1992 ) and 
- much later - to positrons ( Ruiz-Lapuente fc Spruit 19981 : 
iMilne et al.lll999h . as already mentioned above. This leads 
to a decline of the light curve steeper than expec ted for a 
scenario with a constant energy-deposition fraction. I AxelrodI 
( 1980) was the first to describe the 'infrared catastrophe', 
which occurs when the temperature in the ejecta drops be- 
low a critical value. Thereafter, the cooling is mostly accom- 
plished by fine-structure lines in the mid and far infrared. 
While this is merely a re-distribution of the emission and 
has no effect on the true bolometric luminosity, it would 
show up as a faster decline of the UVOIR light curves con- 
structed from observations with a limited wavelength range. 
Dust formation within the ejecta - though a different phys- 
ical process - acts in a similar manner. Conversely, a slow- 
down o f the bolometric decline ra te is observed if 'freeze-out' 
occurs (IFransson fc Kozma|[l993l ). i.e., if additional heating 
is provided by the delayed recombination of electrons at late 
phases. Interaction of the ejecta with circumstellar material, 
resulting in a transformation of kinetic energy into light, 
can have a similar effect. The same is true for possible light 
echoes, where SN light from an earlier epoch is scattered off 
inter s tellar dust clouds towards the observe r (jSchmidt et al.l 
ll994l : ISparks et al.lll999l : iPatat et al.ll2006l ). 

To check for the presence or absence of these effects, 
it is important to model the contribution of the various ra- 
dioactive decays to the bolometric luminosity evolution as 
accurately as possible, including the hitherto neglected IC 
and Auger electrons whose contributions become important 
at very late epochs. To illustrate the effects of these elec- 
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Table 1. Radioactive decay energies (keV decay ^) 



Nucleus 


Auger e 


IC e- 


e+ 


X-ray 


s^Co 


7.594 


10.22 


0.000 


3.598 


56 Co 


3.355 


0.374 


115.7 


1.588 


55 pg 


3.973 


0.000 


0.000 


1.635 


44 Ti 


3.519 


7.064 


0.000 


0.768 


44Sc 


0.163 


0.074 


595.8 


0.030 



trons, we consider the following four decay chains. 



'Ni 



'Ni 



'Co 



*Ti 



*l/2 = 


6.08d 


*l/2 = 


35.60h 


'l/2 = 


17.53/1 


*l/2 = 


58. 9y 



'Co 



'Co 



'Fe 



*Sc 



56 



'l/2 = 


271. 79d 


tl/2 = 


999. 67d 


*l/2 = 


3.97h 



Fe (1) 

^^Fe (2) 

"''Mn (3) 

^■'Ca (4) 



To circumvent the more complicated treatment of the 
photon transport, in the subsequent discussion pertinent to 
very late-time bolometric light curves of different classes of 
supernovae, we limit ourselves to comparing leptonic (and 
X-ray) energy injection rates. We do not include the energy 
produced by the pair annihilation and further assume that 
the kinetic energy of the leptons is completely thermalized 
in situ. Thus, possible escape of positrons (or fast electrons) 
is neglected. 

The energy generation rates presented here do not in- 
clude any heating due to 7-rays and therefore do not accu- 
rately predict bolometric light curves. Our approach how- 
ever does allow for a direct comparison of the relative con- 
tributions to the heating of the positrons produced in the 
decays of ^^Co and **Sc and the electrons produced in the 
decays of ^^Co and ^^Pe. The relevant energies of the dif- 
ferent decay channels are listed in Table [T] These data were 
extracted from the Chart of Nuclides database. National Nu- 
clear Data Centei0. 

3.1 SNe la 

To date, the thermonuclear SNe with data cover- 
age beyond 500 days ar e SNe 1991bg llTura.tto et al.l 
1 19961). 1992 A and 1999by (jCappellaro ct al. 1997*). 2000E 



l Lairetal.1 12009). 2000cx tSoUcrman ct al. 2004) and 
2003hv (iLeloudas et al] l2009l)PrOnly three - SNe 2000cx, 
2001el l|Stritzinger fc SollermanI |2007| ) and 2003hv - have 
been observed in the near-IR at very late epochs. 

SNe la lack an extended envelope, and thus become 
transparent to 7-rays relatively early. In fact, they enter 
the positron- dominated phase about 150-300 days after 
the explosion (iLeiburidgut fc Pintolll992l : iMilne et al.ll200ll : 
ISoUerman et al.l 2004). Between 300 and 600 days, their 



UVOIR light curves seem to follow the Co decay (e.g. 
IStritzinger fc Sollermar] |2007| ). indicating that the major- 
ity of the positrons are trapped. During this phase, flux is 



^ http://www.nndc.bnl.gov/chart/ 

^ For two other SNe la with observatio ns extending to 
~2000 days, SN 1991T l lSchmidt et al.lll994) and SN 1998bu 
llCappellaro et all 120011 ). light echoes dominate the late-time 
emission after ~ 500 days. 




500 600 700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000 
lime [days] 

Figure 1. Instantanous energy gen eration rates fo r initial abun- 
dances taken from the W7 model of lw amoto et all ([1999). Thick 
lines are due to positrons, IC and Auger electrons alone. The thin 
lines also include the full X-ray dose, which is small even in the 
limit of complete X-ray trapping. The contribution of electrons 
and positrons from the decays of 44 -pj g^^j 44 small to be 

seen in this figure. 



progressively shifted from the optical to the near-IR bands, 
the JHK light curves being flat or declining very slowly. 
Starting at 600 days there are some indications of a slow- 
down in the V band below the rate expected for ^^Co decay 
with f ull positron trapping (ISollerman et al.ll2004l : iLair et al] 
l2006l : iLeloudas et alj fioogl ). It has been suggested that a 
heating source other than ^^Co, possibly other radioactive 
species or electron recombination, may be needed to explain 
this slow-down. We point out that a slow-down in the light 
curve is expected when the leptonic energy injection from 
the decay of ^''Co becomes a significant contribution. 

To illustrate this, we plot the leptonic energy generation 
rates of important long lived isotopes between 500 and 2000 
days for th eoretical predictions o f the yields of the fiducial 
W7-model jlwamoto et al.l Il999l ) for a Type la supernova 
(see Fig. [T]). For this particular choice of yields, the slope 
of the light curve begins to deviate appreciably from pure 
56 Co decay after about 800 days. After about 1000 days, 
the heating due to the electrons produced in the decay of 
^^Co equals the heating from the decay of ^®Co (mainly due 
to the positrons). It is also noteworthy that the decay of 
55Fe contributes significantly after about 1400 days (for this 
model), but remains below the contribution of ^'^Co until 
about 2000 days. The contribution of "^"^Ti and its daugh- 
ter '*''Sc is below 10^^ erg s~^ and remains subdominant for 
many years. In reality, the time at which these effects man- 
ifest would depend on the isotopic ratios synthesized in a 
particular explosion. Nevertheless, we predict that a slow- 
down must eventually occur for reasonable nucleosynthetic 
yields. 

Electron capture radioactive decay also involves the 
emission of characteristic fluorescence X-rays (for energies 
see Table [1} ■ Generally, these X-rays are also neglected as 
an energy source for the light curve. Since the opacity of 
the ejecta to X-rays is much higher compared to 7-rays, 
some part of this X-ray radiation may be trapped; there- 
fore this contributes to the heating. To properly quantify 
this would require detailed radiation transport simulations 
including realistic X-ray opacities. However, even under the 



4 Seitenzahl et al. 2009 



extreme assumption that the X-rays deposit all their energy 
(see thin lines in Fig. their contribution is small and our 
conclusions are not strongly affected. 



3.2 SN 1998bw 

SN 1998bw is the first known exampl e of a SN associated 
with a 7-ray burst (|Galama et al.|[T998l) . Spectroscopically it 
lacked any evidence of H or He, and therefore it was classified 
as a SN Ic. It proved to be an unusually energe tic explosion 
of a r nassive, stripped stellar core 10 Ma-,: |Patat et ahl 
I2OOII : iNakamura et al.l l2001bl : iMaeda et all I2OO2I ). Strong 
asy mmetry was inferre d from an analysis of nebular spec- 
tra (|Maeda et al.l [200^ 1. This asymmetry together with the 
high ejecta velocities and the lack of an extended envelope 
help to reduce the 7-ray opacity, so that by 1000 days al- 
most al l 7-rays escape freely (the models of lNakamura et al] 
l|2001al 'l suggest that 1998bw enters the positron dominated 
phase after ~ 400 days). 

Core collapse supernovae synthesize isotopes in differ- 
ent ratios c ompared to SNe la, in p articular they produce 
more ''''Ti. INakamura et al.l (|2001lJ ) calculated nucleosyn- 
thetic yields for hypernovae like 1998bw. Their 30 Bethe 
explosion energy, 10 Mq He-core model predicts a ^^Ni to 
'^''Ni mass ratio 7^57/56 ^ 0.0366, ~ 1.5 time s the solar valu e 
for the '^■^Fe to ^®Fe ratio (which is 0.0234; lLoddersll2003l) 
The is otopic ratios for the asymmetric model of lMaeda et al.l 
l|2002h are similar, albeit wi th a slightly higher contribution 
of ^'*Ti.|H6lich et al.l 11999') proposed an alternative, highly 
asymmetric explosion scenario that requires only 0.2 A/0 of 

After modeling the UVQIR Ught curve to ~ 1000 days 
after the explosion, ISollerman et al.l (I2OO2I ) showed that a 
simple model without contributions from freeze-out effects, 
circumstellar interaction, accretion by a central compact ob- 
ject or light echoes requires a ^^Ni to ^®Ni ratio ~ 13.5 times 
greater than solar. This implies a discrepancy of almost one 
order of magnitude between the value of T^^''/^^ expected 
from explosive nucleosynthesis and that suggested by the 
light curve modeling. However, as we now demonstrate, this 
discrepancy reduces if IC and Auger electrons are included 
in the light-curve calculations. 

In Fig. [2] we show the combined leptonic and X-ray lu- 
minosity c orresponding to y ield pr edictions for the afore- 
mentioned INakamura et al.l l|2001bl ) model. The non-7-ray 
heating due to ^^Co is the dominant contribution between 
~ 1000 and 1600 days. It is remarkable that this very sim- 
plistic approach naturally reproduces the observed slow- 
down of the light curve at ~ 900 days (see Fig. (2)1 without 
the need to invoke any abundance enhancements. 



4 DISCUSSION AND CONCLUSIONS 

We conclude that the electrons produced in the decay of 
^^Co constitute a significant energy source for both stripped 
core collapse and thermonuclear supernovae and should be 
included in codes that model their late-time light curves. 
Fortunately, the vast majority of the energy injected via 
these processes is in the form of primary IC / Auger electrons 
with energy > IkeV. Since the energy deposition channels 
from such high energy electrons are independent of the input 
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400 600 800 1000 1200 1400 1600 1800 2000 
lime [days] 

Figure 2. Instantaneous energy generation rates due to ra- 
dioactivity for initial abundance s taken from the 30 Be thc 10 
Mq He-core hypernova model of lNalcamura et al.l l l2001bl ). Only 
leptons and X-rays (which are small compared to the leptonic 
contribution) are included (i.e. no 7-ray contribution or freeze- 
out effects are conside red). Squares are the data of 1998bw of 
ISoUerman et al] l l2002h . Arrows represent Scr upper limits. There 
are no free parameters. 



electron spectrum (|Kozma fc FranssonI Il992l ). this means 
that standard methods for treating e.g. primary Compton 
electrons are directly applicable. 

We suggest that owing to the neglect of the heating due 
to electrons described in this paper, masses previously de- 
rived for 4*Ti and ^^Ni from light curves should be regarded 
as upper limits. In particular, a reasonable production factor 
-^57/56 jiaturally explains the slow-down in the lightcurve of 
1998bw without the need to invoke extreme super solar val- 
ues. In addition, we predict that after ~ 1400 days, the decay 
of ^""Fe causes a further slow-down of SNe la light curves, 
even though it proceeds via a ground-state to ground-state 
transition. In a scenario with a large escape fraction of the 
more energetic positrons, the effect may be noticeable even 
earlier. 

We limit our discussion in Section 3 to particular classes 
of SNe in which the ejecta are expected to be nearly optically 
thin to 7-rays at the relevant epochs (~ 1000 days) so that 
leptonic heating rates are likely to dominate. However, ^'^Co 
electrons may also be relevant to more complex cases, in- 
cluding SNe II. SN 1987A (Type H-P) is the only SN with a 
bolometric light curve extendi ng from the earliest stages all 
the way to the remnant phase jSuntzeff et alj[l992l : ISuntzefj 
199711). The isotopic ratios of the important radionuclides for 
SN 1987A are not very different fro m the predictions for a 
hypernova (|Nakamura et al]|2001lj ) and the differential ef- 
fect of neglecting the leptonic energy deposition from ^^Ni 
must act in the same sense (i.e. it will lead to an overestima- 
tion of Ti?''^"'^). However, this effect is complicated by the 
large envelope mass, which means the eject a are not fully 
transparent to 7-rays at the relevant epoch jWooslev et alj 
Il989f ) and thus leptonic contributions play a smaller role in 
heating the ejecta. To quantify this requires detailed calcu- 
lations involving both 7-r ay transport and freeze-out effects 
(jFransson fc Kozma|[l993l ). which we defer to future work. 

The additional heating channel pointed out in this 
paper allows, in principle, for a larger escape fraction of 
positrons than previously suggested by modeling late-time 
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light curves - as noted in Section 1. Past analyses of light 
curves have deduced that the escape fraction must be very 
small, but any additional heating source relaxes this con- 
straint somewhat. This may have implications for under- 
standing the origin of the Galactic population of positrons 
and the strength of the Galactic 511 keV annihilation line. 

Our calculations have used atomic data for neutral 
atoms. Although matter in supernova remnants is ionized, 
for light or moderate ionization states, we do not expect the 
IC and Auger decay energies to change significantly, since 
most of the physics concerns the K, L, and inner M shell. For 
very high levels of ionization, the numbers listed in Table [T] 
are inappropriate. However , high levels of ioniza tion are not 
expected at late times fe.g. lKozma fc Franssonlll99& ) unless 
the ejecta is shock-ionized by interaction with circumstellar 
material, a scenario which would, in any case, likely pro- 
hibit reliable abundance determination because of the need 
to separate the contributions from the shock heating and 
radioactive decays to the bolometric light curve. 

The case of SN 1987A (not a very luminous object, 
but at a distance of only ~ 50 kpc in the LMC) demon- 
strates that SN explosions of any type exploding within the 
Local Group can be followed for many years after the ex- 
plosion. Adopting an average rate of ^0.9 SNe (lOOyr)" ^ 
(10^° in spiral galaxies (|Cappellaro et all 1 19991 ) . 

and a Local-Group i?-band luminosity of ^ 6 x 10^° 
l|de Vaucouleurs et al.l Il976l : ISandage fc TammannI Il98ll ). 
statistically a few Local-Group SNe are expected every 100 
years. It is therefore plausible that the effects discussed in 
this paper can be observed unambiguously in the future. 
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